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Abstract. The effect of curvature on the results of frac- 
tal analyses of the galaxy distribution is investigated. We 
show that, if the universe satisfies the criteria of a wide 
class of parabolic homogeneous models, the observers mea- 
suring the fractal index with the integrated conditional 
density procedure may use the Hubble formula, without 
having to allow for curvature, out to distances of 600 Mpc, 
and possibly far beyond. This contradicts a previous claim 
by Ribeiro (1995) that, in the Einstein-de Sitter case, rel- 
ativistic corrections should be taken into account at much 
smaller scales. We state for the class of cosmological mod- 
els under study, and give grounds for conjecture for others, 
that the averaging procedure has a smoothing effect and 
that, therefore, the redshift-distance relation provides an 
upper limit to the relativistic corrections involved in such 
analyses. 
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1. Introduction 

Standard cosmology is based on the assumption that the 
universe is spatially homogeneous, at least on scales suffi- 
ciently large to justify its approximation by a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model. The high 
isotropy measured in the cosmic microwave background 
radiation (CMBR) is usually taken as strong evidence 
in support of this hypothesis. Coupled to a global, re- 
spectively local, Copernican principle, this observation 
has led to the formulation of the Ehlers-Geren-Sachs 
(Ehlers, Geren, Sachs, 1968), respectively almost Ehlers- 
Geren-Sachs (Stoeger, Maartens, Ellis, 1995), theorems. 
Although the assumptions underlying these theorems have 
been discussed by some authors (Clarkson & Barrett, 
1999; Celerier, 2000a), they are generally considered as 
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robust support for the Cosmological Principle. 

The structures seen in galaxy catalogues - groups, clus- 
ters and superclusters, distributed along voids, filaments 
and walls - are not viewed as contradicting this Princi- 
ple, as the scales on which the universe is assumed to 
be homogeneous are much larger than those subtended 
by these structures. Furthermore, the most probable ex- 
istence of dark matter implies that the luminous matter 
seen in the galaxies does not represent the entire source of 
gravitational energy. Thus, the high isotropy of the CMBR 
could, in principle, be reconciled with an inhomogeneous 
luminous matter distribution, provided the missing mass 
for an homogeneous pattern would be provided by dark 
matter. This important issue, which has been formulated 
as the problem of finding whether "light traces mass" is 
still the subject of a number of interesting works (see e.g., 
Bahcall et ah, 2000), but appears far from being settled. 
However, a majority of cosmologists seems to feel more 
confortable with the idea that the structure distribution 
should present a transition towards homogeneity at a more 
or less large scale. 

However, the consensus on a homogeneous feature 
of structures, even on very large scales, has never 
been complete. The debate extends back to New- 
ton (pro-homogeneous) versus Kant and Lambert (pro- 
hierarchical) . The pro-hierarchical approach was reconsid- 
ered in a more modern way by Charlier (1908, 1922), and 
later developed by de Vaucouleurs (1953, 1970), who ex- 
hibited a universal density radius power law. Recent ideas 
on fractals (Mandelbrot, 1982) led to the pioneering work 
by Pietronero (1987). Since then, improvements in obser- 
vational techniques have increased the number and depth 
of available three dimensional galaxy catalogues. A pro- 
gram of systematic analyses of these catalogues has led to 
the claim that one can identify a fractal distribution up to 
150 hr 1 Mpc (Sylos Labini, Montuori, Pietronero, 1998). 
Although these results are disputed by many authors (see 
e.g., Kerscher, 1999), it is now believed by many others 
that one has to increase substentially the 5 ft." 1 Mpc lower 
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boundry generally accepted in the early eighties for the 
scale of the homogeneity transition (Davis, Peebles, 1983). 
The ongoing sharp controversy about this issue awaits a 
new generation of galaxy catalogues in order to be settled 
(Murante et a!., 1998; Martinez, 1999). 

One of the assumptions retained by almost every au- 
thor in the field is that Euclidean geometry may be used as 
a reasonable approximation, up to the distances reached 
by the current galaxy catalogues (Peebles, 1980; Sylos 
Labini et al., 1998). This statement has, however, been 
disputed by Ribciro (1995), who argued that the usual 
practice of disregarding relativistic effects in the study 
of galaxy distribution could cease to be valid beyond 
much smaller scales than commonly believed. Using the 
Einstein-de Sitter model as a theoretical framework, he 
claimed that this approximation could be misleading for 
the scales probed by the most recent catalogues. 



effect of the averaging procedure mainly used for the com- 
pletion of fractal analyses on the relativistic corrections, 
for peculiar (parabolic) FLRW models. We show that, in 
these models, this procedure has actually a smoothing ef- 
fect and we provide grounds to conjecture that this effect 
could apply to any other kind of (more physical) cosmo- 
logical models. 

This article is organised as follows. In Sect. 2, we give a 
short description of the statistical methods used to quan- 
tify the galaxy distribution. Sect. 3 is devoted to an ana- 
lytical study of the Einstein-de Sitter case and to a discus- 
sion of Ribeiro's results. In Sect. 4, we describe the virtual 
catalog method employed to complete our numerical sim- 
ulations, which are presented, with their results, in Sect. 5. 
A discussion and the conclusion appear in Sect. 6. The an- 
alytical expressions used for the calculations are given in 
the Appendix. 



The role of curvature in a universe exhibiting fractal 
structures has also been explored by Humphrey, Matravers 
and Maartens (1998), using, as a toy model, an inhomo- 
geneous isotropic bubble matching a FLRW background. 
Their conclusion was that Einstein's equations and the 
matching conditions imply either a nonlinear Hubble law 
or a very low large-scale density, provided no selection, 
evolution or cosmological constant effect should come into 
play. 

Since then, there has been a growth in the literature on 
this issue. The scale beyond which the relativistic correc- 
tions become non-negligible for the power spectrum anal- 
yses of redshift surveys has, for instance, been recently 
explored by de Laix and Starkman (1998) and Mansouri 
and Rahvar (1999). These authors also claim that a light- 
cone effect is noticeable at lower redshifts than usually 
believed. 

The purpose of the work we present here is to delve 
deeper into the tests of the sensitivity to relativistic cor- 
rections of the methods used in calculating the distribu- 
tion of galaxies. We show that Ribeiro's claim, applied 
to the Einstein-de Sitter case, is ill-founded. Extending 
our study to other homogeneous cosmological models, we 
show that the Euclidean approximation can be considered 
as reasonable up to scales larger than currently reached, 
for a wide class of standard universes, and that, for this 
peculiar class, the redshift-distance relation applying to a 
given cosmological model provides an upper limit to the 
relativistic corrections involved. 

In fact, we know that the Friedmann relation between 
redshift and comoving distance has little chance of be- 
ing valid in our clumpy universe (see e.g., Kurki-Suonio & 
Liang, 1992; Mustapha et al., 1998; Celerier, 2000b). How- 
ever, in this work, we choose, for simplicity, to study the 



2. Quantifying the galaxy distribution 

In this section we give a short description of the methods 
used for quantifying the galaxy distribution. The purpose 
of this digression is to make the arguments concerning our 
main theme clearer. 

The standard procedure, pioneered by Peebles (1980), 
was to follow the methods used for gases and liquids. The 
one-body density, n, is assumed to be constant. The two- 
body density, denoted by n^(r' ,r"), is given by 

n^{r',r") =7i 2 5 (||r"-r'||) 



(1) 



Here <?(r) is the pair distribution function. The main 
interest lies in the deviation of from independent par- 
ticle behaviour and so a two point correlation function 
is defined by £(r) = g(f) — 1. According to the results 
obtained by Peebles, 



ro 



r < r . 



(2) 



The value for ro was given as 5 ft. 1 Mpc. 



It was pointed out by Pietronero (1987) that the ro 
value is a spurious result due to the assumption that n 
is a constant. From Eq. ([!]), we see that the density at r" 
under the condition that there is a galaxy at r' is given by 
the conditional density (r")g(\\r" — r'||). If all galaxies 
are equivalent, we can choose r' — and then, 



>» (1) (r")j(r") = 



i.( 2 ) 



(0,r") 



iW(0) 



(3) 



Assuming isotropy (i.e. all directions to be equivalent), 
the right hand side of Eq.(|J) depends solely on ||r"||, and 
we can write: 



n«(r) ff (r)=r(r), 



(4) 
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Here T is denned by the right hand side of Eq.(||). 
As it has been shown by Pietronero (1987), T(r) is the 
physically significant expression, rather than £(r). If we 
can assume — const., then we can write a relation 
between £(r) and 



e(r) 



r(r) 



1. 



nW{0) 

Just as in Eq.(||), we can write for T, 
r(r) = const.r~~~ tT . 



(5) 



(6) 



It was pointed out by Thieberger et al. (1990) that the 
analysis of data does not necessarily give the same value 
for 7 and 7r- Taking a set of points from a known fractal 
and performing these calculations showed that the theo- 
retical value corresponds to that obtained from yp. 

To reduce the problem caused by finite data sets, 
a method of coarse graining has been introduced, with 
the integrated conditional density T* (Coleman and 
Pietronero, 1992): 



r* = i / Tdv. 

v Jv 



(7) 



It produces an artificial smoothing of rapidly varying 
fluctuations, but correctly reproduces global properties. 
The integral over T is strongly related to the correlation 
integral (Grassberger and Procaccia, 1983), 



C 2 (r) 



1 

iV 7 



1 



N 



(8) 



O is the Heaviside function. The inner summation is 
over the whole set of N — 1 galaxies with coordinates X.,- , 
j i, and the outer summation is over a subset of N' 
galaxies, taken as centers, with coordinates Xj. By taking 
only the inner N' galaxies as centers we allow for the effect 
of the finiteness of the sample, (see e.g. Sylos Labini et al., 
1998). This procedure is also widely used in astronomy 
(Provenzale, 1991). 

3. Analytical study of the Einstein-de Sitter case 

As Ribeiro (1992, 1995; hereafter refered to as R92 and 
R95) and other authors cited therein quite rightly pointed 
out, galaxy observations are carried out along our past 
light cone. In a homogeneous FLRW cosmological model, 
the proper density is constant on hypersurfaces of con- 
stant proper time. Therefore, as one looks back to higher 
and higher redshifted galaxies, the crossing of the null 
geodesies through hypersurfaces of constant proper time 
and density implies changes in the observed density and an 
apparent inhomogeneity can be measured. One can thus 
wonder whether the inhomogeneous pattern identified by 
the fractal analysis of galaxy counts is a mere artefact of 



this light cone effect. 

In R95, the magnitude of this effect in an Einstein-de 
Sitter universe is examined. The conclusion is that "even 
accepting an error margin of 25% in the measurements 
of the global density, a redshift equal to 0.1 is approx- 
imately the deepest scale where we could observe a ho- 
mogeneous distribution of dust in an Einstein-de Sitter 
model", unimpaired by the light cone effect. Now, such 
a redshift corresponds to a luminosity distance of about 
410 Mpc for a value of the Hubble constant Ho = 75 Km 
s" 1 Mpc" 1 and 310 Mpc for H = 100 Km s" 1 Mpc" 1 . 
This is actually less than the depth reached by the ESO 
Slice Project galaxy survey (Vettolani et al., 1997, 1998), 
which recently fueled the controversy about the transition 
towards homogeneity of the luminous matter distribution 
(Scaramella et al., 1998; Joyce et al., 1999). Other surveys 
aimed at going beyond this limit are currently under way. 

Ribeiro's argument is based on the analysis of the sen- 
sitivity of the integrated conditional density T*, Eq.(^), 
relative to the redshift. We show, in the present section, 
that the expression retained by Ribeiro for T* has to be 
modified, which results in T* being less sensitive to the 
redshift value than claimed in R95, thus impairing its con- 
clusions. 

In R92, the author establishes expressions for some ob- 
servational quantities, in an Einstein-de Sitter space-time, 
as a function of the radial comoving coordinate r, chosen 
as a parameter along the null geodesies. We list below, and 
give in R95 units, the quantities relevant for our purpose. 

The cumulative number count N c (r) is the number of 
sources which lie at radial coordinate distances less than 
r, as seen by the observer at r = 0. In R95, it is written 

as 



N B = 



2r J 



9Mg ' 

where Mq is the average galactic rest mass. 



(9) 



The luminosity distance d\ of a source is the distance 
from which the radiating body, if motionless in an Eu- 
clidean space, would produce an energy flux equal to the 
one measured by the observer. It thus verifies 



/ = 



(10) 



L being the absolute luminosity, i.e. the luminosity in the 
rest frame of the source, and / the measured bolomctric 
flux, i.e. integrated over all frequencies by the observer. 
Its R95 formulation is 



d, =9r 



3-ffn 



18 



(11) 
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And the redshift z runs as 



, 18 

1 + *='«o 



18 



(12) 



From the definition of T, given in the previous sec- 
tion, one obtains, in Euclidean space-time, (see also Sylos 
Labini et al., 1998): 



r(r) 



1 dN r 



S dr ' 

where S is the area of a spherical shell of radius r. 



(13) 



Arguing that the relevant distance for an analysis per- 
formed in a curved space-time is di, Ribeiro defines r(d;) 
as 

1 dN r 



Sd{di) 



(14) 



where S is presented as being the area of the observed 
spherical shell, of radius di, and is thus written 



S(di) = iirdf 



(15) 



It is here worth noting that the luminosity distance is 
the observable quantity relevant for radially measured dis- 
tances. It is thus correct to use this quantity to evaluate an 
elementary change in the number density N c as measured 
by the observer on his light cone. But, when looking at a 
cross sectional area dS perpendicular to the light ray and 
subtending a solid angle eifi, the observer must consider 
the area distance d a , also called "observer area distance 
" by Ellis (1971) and "corrected luminosity distance" by 
Kristian and Sachs (1966), defined by 



dS = d„dQ, 



(16) 



and related to the luminosity distance of an object mea- 
sured with the same redshift by 



di = d a (l + z) 1 



(17) 



Therefore, the area of the observed spherical shell, at 
a given redshift z, used in Eq.(|l4]) must be written 



s = And * 



{l + zf 



(18) 



From Eqs.(||) and we get an expression for di as 
a function of the redshift z: 



-Wo 



(19) 



T(z) 



Using Eqs.(g), (gj), (yj) and (|19|), we then derive T(z 
3H$ ( l + z 



(20) 



8ttM g \2y/TTz-l, 

As can be easily seen, by comparing Eqs.([l5|) and (|l8|), 
this expression differs from the one that would proceed 



from R95 by a factor (1 + z) 4 . Therefore, its sensitivity to 
the redshift value is much lower. 

The volume element dV, relevant for calculations on 
the light cone, is 



dV = Sd(d t ) = 4tt- 



rd(di 



(i-M— (21) 

Substituting Eq. (|l9|) into Eq.(|2l|) and integrating, we 
obtain 



V = 



32ir (yT + I - l) 3 (2z + 6yT + I - 1) 
H§ 15(1 + z)i 1 j 

Substituting this expression into Eq. (Q) and using 



Eq.(£0D, we get 

T*(z) 



15# n 2 



8irM G (2z + 6vT+T - 1) 



(23) 



which can be easily compared to Eqs.(19) and (20) in R95. 
We thus see that, contrary to Ribeiro's claim, T* differs 
from the average number density < n > if the correct rel- 
ativistic quantities are retained. It is interesting to note 
that, as di(z), T* is an increasing function of z, not a de- 
creasing one, as is Ribeiro's < n > (z) function. 

Furthermore, if we take temporarily Hq = 75 Km s _1 
Mpc -1 to compare our results with those proposed in R95, 
we see that between z = and z = 1 the increase of our T* 
is only 5% (instead of a 87.5% decrease for Ribeiro's). A 
significant departure of T* from its value in an Euclidean 
universe should therefore occur at redshifts higher than 
claimed in R95. Its ensuing conclusions thus do not agree 
with the analysis, which is confirmed by the results of our 
simulations presented below. 

It is also interesting to note that the 5% increase in 
our r* at z = 1 is markedly less than the 17% increase 
of the di/rE ratio, where te = z/Hq (in units c = 1) is 
the Euclidean distance. In an Einstein-de Sitter model of 
the universe, the redshift-distance relation therefore pro- 
vides a large upper-limit for the departure of T* from its 
Euclidean counterpart. The simulations presented below 
show that this result can be generalized to a wide class of 
standard homogeneous models of the universe. 



4. The virtual catalog method 

A real catalog gives, for each observed galaxy, a set of four 
quantities: two angles 9 and <f>, the redshift z (or, equiva- 
lently, the recession velocity v) and the apparent luminos- 
ity / (or, equivalently, the apparent magnitude m). This is 
usually obtained by measuring galaxies with / > fu m (or, 
equivalently, m < mum), in a region of the sky delimited 
by a solid angle f2. 



M.N. Celerier & R. Thieberger:Relativistic corrections to fractal analyses of galaxies 



5 



In an analysis as described in Sect. 2, two of the above 
observablcs arc transformed, for each galaxy, into non- 
observable quantities, the Euclidean distance te and the 
absolute luminosity L (or, cquivalently, the absolute mag- 
nitude M), using the following equations: 



cz 
So' 



L = A-Kr 2 E f, 

or, equivalently, in Mpc: 
M = m — BlogiorE — 25. 



(24) 



(25) 



(26) 



Two different kinds of galaxy catalogues are currently 
available. Pencil beams are generated by deep surveys of 
a region of the sky delimited by a very narrow solid an- 
gle fl. Much wider distributed samples are generally more 
limited in depth. A peculiar analysis method applies to 
each of these two types of datasets. 

Owing to their very small orthogonal extension, pen- 
cil beam catalogues can only be analysed using the ra- 
dial count iV c (< r) from the observer (Joyce et al., 1999). 
As we stressed in Sect. 3, the sensitivity of the cumulative 
number count iV c to relativistic corrections is well rep- 
resented by the redshift-luminosity distance relation. In 
the Einstein-de Sitter model, the departure of this rela- 
tion from its Euclidean counterpart is only 17% at z = 1, 
which is less than the 25% error margin in the measure- 
ments accepted in R95. The Euclidean approximation can 
thus be considered as valid up to higher rcdshifts than 
z ~ 0.3, currently reached in these surveys. In the follow- 
ing, we shall not be concerned with the study of relativistc 
corrections applying to this method, but shall focus our 
attention on the more accurate one described below. 

The second method, in which one averages over a large 
number of (inner) galaxies, allows us to fit, inside the sam- 
ple volume, spheres with a sufficiently large maximum ra- 
dius, since the T statistic can only be computed up to 
scales limited by this radius. To avoid selection effects, one 
generally extracts from the data a set of volume-limited 
(VL) samples, each sample containing every galaxy with 
an absolute luminosity below Lvl (or, equivalently, with 
an absolute magnitude beyond Mvl) and distance below 
Rvl, such as 



Lvl = ^R VL f 

or, equivalently, in Mpc: 



M VL 



blogioRvL — 25. 



(27) 



(28) 



The fractal analysis is thus applied to the set of VL 
samples (Coleman, Pietronero, 1992; Sylos Labini et al., 
1998). Therefore, relativistic corrections must be consid- 
ered twice. First, for the determination of distances, and 



second, for the construction of the volume- limited samples 
extracted from the catalog. As expected, our numerical 
simulations show that the dominating effect is the dis- 
tance (or light cone) effect. 

We here propose to approach the study of these rela- 
tivistic corrections from a different point of view, by con- 
sidering the light cone problem stated as follows. 

A homogeneous universe, with a random galaxy dis- 
tribution on each constant proper time hypersurface, ex- 
hibits, when observed on our past light cone, a given 
pattern of structures (see above Sect. 3). If, for a chosen 
FLRW model of the universe, we can construct, from a 
random distribution of points, a virtual catalog as it could 
be actually observed, we can apply, to this virtual catalog, 
the above described classical statistical analysis method, 
in the Euclidean approximation, and consider the fractal 
dimension obtained. The departure of this dimension from 
its homogeneous value, D — 3, provides a measure of the 
error made when ignoring curvature effects. 

In this work, we limit ourselves to the study of spa- 
tially flat FLRW models. Spatial flatness seems to be a not 
too restrictive prescription, as it is compatible with the 
most recent results of CMBR anisotropy measurements 
(de Bcrnardis et al., 2000; Balbi et al., 2000). 

The line-element in proper time t and comoving coor- 
dinates r, 9, <p, can thus be written 

ds 2 = c 2 dt 2 - a 2 (t)[dr 2 + r 2 {d0 2 + sin 2 9d<p 2 )] (29) 

where a(t) is the scale factor which, in the matter domi- 
nated area, evolves as 1 3 . 

On each hypersurface t = const., the galaxy distri- 
bution is random (homogeneous), with an absolute lumi- 
nosity distribution given by an assumed luminosity func- 
tion 4>(L). Galaxies evolve on world lines (r = const., 9 = 
const., (j> — const.). Thus, if we retain the assumption 
that, at the distances probed by the current surveys, 
galaxy evolution can be ignored, a given galaxy, with co- 
moving coordinates (r, 0, (f>) and absolute luminosity L, 
will keep the same values for these quantities on each 
t = const, slice. 

In the geometrical optics approximation, the light re- 
ceived from the galaxies travels on null geodesies (see e.g 
Kristian and Sachs, 1966) and thus the objects compiled 
in the catalogues are all located on our past light cone. 
At each measured value of the redshift z corresponds an 
intersection of this light cone with some t = const, hyper- 
surface. 

To reconstruct the catalog which would be observed, 
in a given FLRW universe, by an earth-grounded observer 
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looking on our past light cone from the "here and now" 
point (to, ro = 0) to a random distribution of galaxies, we 
adopt the following method. 

First, on the "now" hypersurface t — to, we generate a 
random distribution of points, i.e. galaxies, fixing, for each 
of them, values for r < r max , 6 < 9 max , 4> < <p max and 
L > L min . Contrary to the producers of real catalogues, 
we are not limited in angular observation, so we can take 
the angular maximum values as large as we wish, up to 
a complete sky. The absolute luminosity distribution is 
given by the luminosity function retained, see Sect. 5. This 
t = to hypersurface corresponds to a spatial section of a 
universe which has evolved from different sections t < to 
up to t = to, to reach this random configuration. 

Then, for each galaxy with a radial comoving coordi- 
nate r and absolute luminosity L, as first generated on 
the t = to hypersurface, we calculate the redshift z and 
apparent luminosity /, as they would be measured on the 
observer's past light cone, if the chosen FLRW model was 
to constitute a reliable description of our universe. For this 
purpose, we use the definition of the luminosity distance: 



/ = 



L 



4ndl 



(30) 



where we substitute the expression of di as a function of 
r relevant in the given FLRW model. The redshift z also 
proceeds from its expression as a function of r in this 
model (see the Appendix). 



We introduce a cut-off at low luminosities (we chose 
L = 0.05) and at high luminosities (we chose L = 10.). 
From (L, x, y, z) we construct, as explained in Sect. 4 and 
in the Appendix, (f,z,(f>,9). For the Einstein-de Sitter 
case, we use Eqs.(^3), ( ]35| ) and J3C|). For the parabolic 
homogeneous models with a non-zero cosmological con- 
stant, we use Eqs. (|40|), (El]) and (|o|). Our choice of the 
parameters are 

Hq = 65 Km s" 1 Mpc" 1 

Vt M = 0.3 

Q A = 0.7. 

Our next stage is to obtain the dimension. We use for 
that purpose Eq. (||), which is essentially like using T* . 
As is mentioned in Sect. 2, we have an inner region at a 
distance larger than some T> max from the border of the 
total sample. We now calculate distance and luminosity 
classically, using the Euclidean formulae: 



te = 



L, 



V 



(32) 



(33) 



Our limited sample requires Le > L min . 

To understand better our results we perform an addi- 
tional type of calculation, considering just two points of 
the inner region, closest to the center and farthest from 
the center. For each point we can check the local density 
in addition to the dimension in that region. 



We thus obtain the virtual catalog corresponding to a 
homogeneous galaxy distribution in this model, as a set of 
quadruplets (f,z,9,(f>), each corresponding to a virtually 
measured galaxy. 

Now, we can perform a classical statistical analysis of 
this virtual catalog, as described in Sect. 2. We can thus 
obtain, varying the depth of the data sample, different val- 
ues for the calculated fractal dimension D. The departure 
of D from its homogeneous value, D = 3, gives the mea- 
sure of the error made when ignoring curvature effects, 
at the scale probed by the depth of the retained sample, 
provided the given cosmological model can be considered 
as a good approximation of the observed universe. 



5. Numerical simulations 

5.1. The method of calculation 

We start with random (x,y,z). Next, we discard all the 
points for which R is larger than R max . Then, we assign 
to each point a "Luminosity" , L, obtained from the dis- 
tribution: 



5.2. Results 

For each of the three cosmological models studied, we per- 
formed calculations for six different sets of random num- 
bers and took their average. By using these different six 
catalogues for each cosmological model, we could obtain 
errors on our results. 



We performed the calculations for different L m 
convenience we used a parameter Lm — const * L r 



For 



P(L) = A*exp(-L)/L. 



(31) 



We then compared our results to those we obtained 
when ignoring the relativistic corrections for the same 
random numbers. We compared the dimensions obtained 
in both cases using a method developed by Benzi et al. 
(1995). The advantage of this method is that part of sys- 
tematic errors, like border problems, are less pronounced. 

We checked the case V max = 200 Mpc. Comparing it 
to the non- relativistic case, no significant difference could 
be detected. 

Checking the differences for the closest and farthest in- 
ner point, we obtained that, for the non-relativistic case, 
33,162 points were inside the sphere, when taking the 
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Table 1. Density results for the three models 



Model# 


Pret 


M1/M2 




D(R m ax) 


1 


100 


1.35 ± 0.02 


3.01 ± 0.04 


2.97 ± 0.03 


2 


100 


1.12 ± 0.01 


3.00 ± 0.03 


3.00 ± 0.01 


3 


100 


1.00 ± 0.01 


3.00 ± 0.02 


3.00 ± 0.01 


1 


25 


1.41 ± 0.02 


2.99 ± 0.03 


2.99 ± 0.02 


2 


25 


1.29 ±0.01 


2.98 ± 0.03 


3.01 ± 0.01 


3 


25 


1.22 ±0.01 


3.00 ± 0.02 


2.97 ± 0.03 


1 


8 


1.46 ± 0.02 


2.99 ± 0.03 


3.01 ± 0.02 


2 


8 


1.42 ± 0.02 


2.97 ± 0.02 


2.99 ± 0.04 


3 


8 


1.40 ± 0.02 


2.99 ± 0.02 


3.01 ± 0.03 



- Model #1 is the Einstein de Sitter Model. 

- Model #2 is the Parabolic Homogenous Model 
0.3, Ov = 0.7. 

- Model #3 is the Parabolic Homogenous Model 
0.01, Q A = 0.99 



with !Jm — 
with O.M = 



closest point (R m in = 550 Mpc). For the farthest point 
{Rmax = 1100 Mpc), we obtained, in the same sphere, 
33,152 points. This was expected as we chose a homoge- 
neous sample. 

The results are presented in Table 1. We denote by Ml 
and M2 the number of points in the two spheres of radius 
T^max, whose centers are R m in and R max respectively, and 
by Pret, the percentage of points retained in each sam- 
ple, i.e. for three values of L m i n . We obtain definitely an 
influence on the average density, depending on the spe- 
cific model. For the three models, even though this value, 
(M1/M2), is quite appreciable, still this does not affect 
the dimension within a T> max — 200 Mpc radius. For all 
the cases, we obtain: D = 3.00 ± 0.03 . 

Our next purpose was to check to what extend the 
dimension diminishes when we go to larger distances. 
We calculated the case where we consider the whole 
sphere around our galaxy and take for the inner region 
a midrange of the radius so that the centers in that region 
are farther than T> max from the borders. The calculation 
was performed just for one case, the Einstein-de Sitter 
model, and L min was chosen so that 42 % of the galaxies 
where retained, on average. We averaged over six sets of 
120,000 points, chosen randomly, and over all the above- 
mentioned points of the central region (about 12,000), and 
checked the decrease in the dimension when averaging over 
the ranges 150 -300 Mpc versus 450-600 Mpc. We obtain 
that the dimension slightly decreases by 1.6±0.2%, which 
is consistent with our other calculations. 

To compare these weak discrepancies with the light 
cone effect on the luminosity-distance relation, we note 
that 

- for the Einstein-de Sitter case, at T> max = 200 Mpc, the 
contribution of the second order term in the power series 
expansion of di is a little more than 1% of the total in 



Eq.(|6|), but at V max = 600 Mpc, it is more than 3%. 

- for the case !1m = 0.3 and S1a = 0.7, at V max = 200 
Mpc, this contribution is about 3.2% of the total in 
Eq.©. 

- for Q M = 0.01 and Q A = 0.99, at V max = 200 Mpc, it is 
more than 4%. 

We thus see that, except for the Einstein-de Sitter case 
at T> max = 200 Mpc, this effect is stronger on di than 
on D. Furthermore, as, at T> max — 200 Mpc, this effect 
on the fractal dimension D seems to be model indepen- 
dent, the error made by ignoring the relativistic correc- 
tions in classical analyses must be smaller, at this scale, 
than the statistical noise of the method. Therefore, we can 
safely conclude that, within the studied class of models, 
the redshift-distance relation can actually be used to pro- 
vide an upper limit to the relativistic corrections needed 
to complete this kind of analysis. 



Furthermore, as Eq.(|17|) provides a general relation be- 
tween di and the area distance d a , the smoothing (1 ± z) 4 
term appears in the calculation of T* whatever cosmolog- 
ical model is retained to represent the observed universe. 
We are thus inclined to suggest that the above property 
could apply to any model and that the averaging proce- 
dure has actually some general smoothing result. 

6. Discussion and conclusion 

In the present article, we have studied the effect of curva- 
ture on the results of fractal analyses of galaxy distribu- 
tion. 

First, we explored analytically the Einstein-de Sitter 
model and found that, contrary to the claim in R95, the 
sensitivity of the integrated conditional density T* to the 
redshift is far less than the sensitivity of the luminosity 
distance. Therefore, in this model, the redshift-distance 
relation can be used to derive an upper limit to the rela- 
tivistic corrections applying to classical fractal analyses. 

Then, we enlarged our study to a numerical simula- 
tion of a wide class of parabolic homogeneous models. To 
resolve every ambiguity which could arise from a mixing 
of classical and general relativistic notions of distances, 
surfaces, volumes and averagings, we calculated the error 
which could be made by an observer, immersed in a FLRW 
universe, and analysing, with Euclidean statistical tools, 
an intrinsic homogeneous pointlike galaxy distribution, as 
measured on his past light cone. We have shown that, in 
these models and up to distances far beyond those probed 
by the current analyses, the curvature effect is negligible 
compared to the measurement uncertainties. For Hq = 65 
km s _1 Mpc -1 , the three tested models exhibit, at the 
scale of 200 Mpc, a discrepancy for the numerically calcu- 
lated fractal dimension of on the order of 1%, and a 2% 
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limit still holds for the Einstein-de Sitter case at scales on 
the order of 450-600 Mpc. 

In fact, contrary to the author's claim in R95, the pro- 
cedure averaging over many points, far from increasing the 
curvature effect, has some smoothing results, and there- 
fore this effect is weaker on the fractal dimension than 
on e.g. luminosity distances. We have actually proved this 
property for the class of universes under study and have 
given grounds to suppose that it can apply to any other 
cosmological model. 

Therfore, as regards the current and forthcomimg re- 
sults which could be derived from Euclidean statistical 
analyses of galaxy data, the relativistic corrections could 
be consistently ignored, up to the abovementioned scales, 
and probably far beyond, in every universe for which the 
dependence of the luminosity distance on the redshift does 
not too widely depart from the Hubble law at the probed 
scales. 

We want to mention, for completeness, that the cur- 
vature correction is not the only bias that can impair the 
results of such analyses. Evolution effects and, mainly, K- 
corrections are known to generate other kinds of errors, 
but are beyond the scope of the present work. For a dis- 
cussion of the impact of K-corrections on such an issue, 
we refer the interested reader to Scaramella et al. (1998) 
and Joyce et al. (1999). 
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Appendix A: The FLRW models tested 

We give, in this appendix, the expressions for di and z 
retained for the conversion of the (L, r, 0, 4>) data into the 
(/, z, 8, (j)) quadruplets, in the FLRW models chosen to 
perform the numerical simulations described in Sect. 5. We 
have limited these simulations to three examples of the 
spatially flat subclass: the Einstein-de Sitter case, with 
Qm = 1 and f^A = 0, currently the most popular model 
with — 0.3 and fl\ = 0.7 and a low matter density 
case, with VLm — 0.01 and £1a = 0.99. As it can be easily 
seen from a comparison of the di expansions in powers of 
z calculated below, the low matter density case exhibits 
the larger second order term correction for "realistic" flat 
models with < CIm < 1. It can thus been considered as 
a limiting case for the study presented here. 



For the Einstein-de Sitter case, we use Ribeiro's rela- 
tions in R95, corrected for consistency from a dimensional 
point of view: 



di = r I 1 — — -r 
1 2c 



1 i'-f- 



The power series expansion of di thus reads 



d, = — ^ — -l '"><-'' 



H 4H Q 



0{z i 



(34) 
(35) 

(36) 



For the parabolic homogeneous model with a non-zero 
cosmological constant, the following expressions can be 
found in text books (see e.g. Carroll, Press, Turner, 1992): 



di = (1 + z)r, 



(37) 



d, = 



c(l + z) 
H 



[(1 + z') 2 {l + fl M z') - z'(2 + z')n A ]-2dz'. (38) 



It is easy to derive an exact expression for r from 
Eqs.@ and (§§): 

r =TT I [(lW) 2 (l + fW)-^(2 + z')f> A r^z',(39) 
n o Jo 

For small z, we can Taylor expand the right-hand side 
of Eq. (139) and invert to obtain 



Z = ( V 1 + V (2 " A - °" - 2)r - 1 ) - (40) 

Substituting Eq. jic| ) into Eq. ( |37| ) , we get di as a func- 
tion of r: 



di = 



1 



TT 

1 + — (2fi A - Q M - 2)r - 1 
c 



2J7a — "M ~ 2 
and can derive the power series expansion of di : 

.2 , (42) 



■(41) 



cz c 

di = — — h ——(2 — Qm + 2Q A )z 2 + 0{z 6 ). 
tlo 47i 
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